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TWO NEW ZETA CONSTANTS: FRACTAL STRING, CONTINUED 


FRACTION, AND HYPERGEOMETRIC ASPECTS OF THE RIEMANN ZETA 


FUNCTION 


STEPHEN CROWLEY 


ABSTRACT. The Riemann zeta function at integer arguments can be written as an infinite sum of 
certain hypergeometric functions and more generally the same can be done with polylogarithms, 
for which several zeta functions are a special case. An analytic continuation formula for these 
hypergeometric functions exists and is used to derive some infinite sums which allow the zeta 
function at integer arguments n to be written as a weighted infinite sum of hypergeometric 
functions at n — 1. The form might be considered to be a shift operator for the Riemann zeta 
function which leads to the curious values ¢F (0) = Io(2) — 1 and ¢F (1) = Ei(1) — y which 
involve a Bessel function of the first kind and an exponential integral respectively and differ from 
the values ¢ (0) = -4 and ¢ (1) = oo given by the usual method of continuation. Interpreting 
these “hypergeometrically continued” values of the zeta constants in terms of reciprocal common 
factor probability we have ÇF (0)~! = 78.15% and ÇF (1)! = 75.88% which contrasts with the 
standard known values for sensible cases like ¢ (2)~1 = Ss = 60.79% and Ç (3)! = 83.19%. 
The combinatorial definitions of the Stirling numbers of the second kind, and the 2-restricted 
Stirling numbers of the second kind are recalled because they appear in the differential equation 
satisfied by the hypergeometric representation of the polylogarithm. The notion of fractal strings 
is related to the (chaotic) Gauss map of the unit interval which arises in the study of continued 
fractions, and another chaotic map is also introduced called the “Harmonic sawtooth” whose 
Mellin transform is the (appropritately scaled) Riemann zeta function. These maps are within 
the family of what might be called “deterministic chaos”. Some number theoretic definitions are 
also recalled. 
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1. THE ZETA FUNCTION 


1.1. Riemann’s ¢(t) Function. Riemann’s zeta function, named after Bernhard Riemann(1826- 
1866), is defined by 


C(t) = paan Vite C: R(t) > 1} 
(1.1) = ps eo (2n4+1)* vte Cc: Rit) > 1} 
= (1—27*) n (t) Vite C: R(t) > 0} 


where R(t) and Z(t) denote real and imaginary parts of t respectively and 7 (t) is the Dirichlet eta 
function, also known as the alternating zeta function, named after Johann Dirichlet(1805-1859) 








(1.2) n(t) = Ere Se V{t € C: Rit) > 0} 
= ph le eade V{[te OC: Rif) > 0} 


where the integral is a Mellin transform of (e* + 1)~*. The function Ç(t) is analytic and uniformly 
convergent when R(t) > 1 or R(t) > 0 when using the eta function form. The only singularity of 
¢(t) is at t = 1 where it becomes the divergent harmonic series. The reflection functional equation 
[48, 13.151] which relates ¢(t) to ¢(1 — t) is given by 


(1.3) ¢(t)n—*2' T(t) cos (5) =Ļ(1-t) 


The interpretation of zeta in terms of frequentist probability is that given n integers chosen at 
random, the probability that no common factor will divide them all is ¢(n)~'. In other words, 
given an array i of n random intgers, ¢(n)~‘ is the probabability that gcd (i4, i2,...,in) = 1 where 
gcdis the greatest common denominator function. So for example, the probability that a pair of 
randomly chosen integers is coprime is ¢ Qi" = $ = 60.79%, and the probability that a triplet 


2 
of randomly chosen integers is relatively prime is ¢ (3)~* = 83.19%. [37][48, 13.1][7, 1.4] 


1.1.1. The Generalized Hurwitz Zeta Function ¢(t,a). A more general function which includes Rie- 
mann’s Zeta function was defined by A. Hurwitz. 

co 
(1.4) C(t,a) = X (nta) 

n=0 
Notice that the summation starts at n = 0 whereas Riemann’s starts at n = 1. It is apparent that 
¢(t) is a special case of ¢(t,a) where 


(1.5) 6) = ont = CD = ns 
n=1 n=0 


[14][48, 13.11] 
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1.1.2. Hypergeometric Representations of the Lerch Transcendent: ®(z,t,v). The Lerch transcen- 
dent ®(z,t,v) [10, 1.11] is a further generalization of the Hurtwitz zeta function 


(1.6) O(z,t,v) = eo wan 

valid V|z| < 1 or {|z| = 1 : R(t) > 1} which is related to ¢(t,v) and ¢(t) b 
O(1,t,0) =C(t,v) 

(1.7) ®(1,¢,1) =¢() 


@(1,t,1/2) =C(t,1/2) = (2 — 1)c(t) 


When t = 1 the Lerch transcendent reduces to 





(18) en( igu e) 
®(z,1,v) = > 

and when n €E {0,1,2,....}, ®(z,n,v) has the hypergeometric representation [19] 
(1.9) O(z,n,v) = vi Fn ( i am z) 

+ Un 
yielding 
(1.10) C(n,v) =u, Fn ( : eas ) 

+ Un, 

and 


— (_2" 1 T 
e = (2) a al  S 


Tt 
Fy 
=n+1 (5 am ) 


and thus due to (1.1) and (1.7) we have the hypergeometric transformation 
1 IŻ I 

1.12 pat Be, Oh (Cm Gee? a ew ie (a 

(1.12) e ( ii | aan (2 ) 


where the argument absent in pF} is assumed to be 1 and the symbol ¢,, denotes a parameter vector 
of length n where each daa is equal to c (e.g. 53 = [5, 5, 5]). 


(1.11) 


1.1.3. The Hypergeometric Polylogarithm. The polylogarithm, also known as Jonquière’s function, 
is defined V|t| < 1,n € {0,1,2,...} by 


: oo k 
Lin(t) = $k a a 
ln 
=n+1 Fy, ( il ) t 


Tt (Qe 
(1.13) E mia oR Hee 








(2) 
T(k+1 
= trp st ER 


TOP EN TKY 
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The hypergeometric representation (4.2) of Lin (t) =n4i Fn ( oe e): = Lif (t) where 
1---On 


Q,..-An41 = Inji and bi ...bn = 2n is nearly-poised of the first kind [41, 2.1.1] since a, + bı = 
3=...=a,+b, =3. The notation Lif (t) refers specifically the hypergeometric form of Lin (t). 
The derivatives and integrals of Li,,(t) satisfy the recurrence relations 


: Lipa 
(1.14) m rE 
: 1 1 
d n+1 n+1 
gan(t ona Ee) 
jitas Stima 
1.15 1, ile 
l ) i n+1Fn ( 3 i s) ds = n42Fn+1 ( 5 = t) t 


and the reflection functional equation for Lin (1) = ¢(n) is 


(1.16) ERG = gu} 


Li? (t) is seen to be (n — 1)-balanced (4.3) with the trivial calculation 


n n+1 
(1.17) X2- 1=2n-(n+1)=n-1 
k=l k=1 


The usual defintion of Lin (t) requires analytic continuation at t = 1 but this is not necessary because 
the hypergeometric function converges absolutely on the unit circle when it is at least 1-balanced 
(4.3) which is true Vn > 2. The only Saalschtitzian polylogarithm is Lig(t) [32, Eq3.8] [20, 25:12][26, 
1.4.2] 


1.1.4. The Differential Equation Solved by Liž (t) and Some Combinatorics. Some combinatorial 
functions need to be defined before writing the differential equation solved by Lif (t). Let a partition 
be an arrangement of the set of elements 1,...,& into n subsets where each element is placed into 
exactly one set. The number of partitions of the set 1,...,k into n subsets is given by the Stirling 
numbers of the second kind [2, 1.1.3][38, 2.7] defined by 


k A jk n 
Ge 


1.18 See 
(1.18) = bjo TOFTA- | 

ee 1= n, 2p-1 
= “Tay, fk er 





The ,F;,_1 representation of { : \ is (n—k)-balanced (4.3) since (k—1)—((1—n)+2(k—-1)) = n—k. 


The r-restricted Stirling numbers of the second kind { i \ , or simply the r-Stirling numbers, 


3 
counts the number of partitions of the set 1,...,n into k subsets with the restriction that the numbers 
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1,...,r belong to distinct subsets. [29] The recursion satisfied by { : \ is given by 


s 


0 k<r 
k On,r k=r 

(1:19) vat r Tag fey 
r n + n>r 

n : n—-1 X 


where n,m = { : $ 7 ks is the Kronecker delta. Specifically, the 2-restricted Stirling number- 


s[15, A143494] appearing in the differential equation for Lif (t) are given by 
k O fk k—1 
n f, inf ln 
k- kf k-2 : ie 
(x20) = et D #7? gra 


k-32 (j+2)"7?(-1)? 
(-1)* yo 5-0 aT 


II 


The (n + 1)-th order hypergeometric differential equation (4.5) satisfied by f(t)=Lif (t) (1.13) 








fO+GZfOC - 4) n=0 
1.21 0= n We n+1 ae n+1 
HF) + Dont afer) (erry BT bowel BAT aa 
2 
has a most general solution of the form 
(1.22) F (2) =2 + yGal(t) + EEZ zm Int)” 
where 2, Yy, 21,---,2n—1 are arbitrary parameters and G,(t) satifies the recursion 
ih n=0 
(1.23) Gol) E a S E n=2 
f =at n>3 
which has the explicit solution 
= n= 
In(t p= 1) K a e 
(1.24) Gut) =$ Li(1—t) + In(t — 1) n(t) n=2 
eee rn In(t)"~1 + 7 BO Li, (t) n>2 
The indicial equation of (1.21) at the t = 1 is 
—1)” 71T (n—1— =n 7 
(1.25) ind(Lif (¢)) = -= PO ee 


The (n + 1) roots of ind(Li* (t)) are the exponents of (1.21) which are simply 
(1.26) {t : ind(Lif (t)) = 0} =0,1,...,.n-1,n-1 


where the last root n — 1 of ind(Lif (t)) is the balance of Lif (t) (1.17) having multiplicity 2 thus 
inducing the logarithmic terms of (1.22). [16, 15.31 and 16.33] These equations were derived by 
writing the differential equation for increasing values of n and then noticing that the developing 
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pattern of coefficients were combinatorial. After deriving the general combinatorial differential 
equation, it was solved for increasing values of n which resulted in nested integrals of prior solutions 
and then the general solution was derived from that pattern. 


1.1.5. The “Hypergeometric Form” of the Zeta Function. The main focus will be on the special 
case Lif (t) at unit argument where it coincides with the Riemann Zeta function at the integers. 
As with Li’ (t), the symbol ¢/"(n) refers specifically to the hypergeometric representation of ¢(n) 
at non-negative integer values of n. Using (1.5) and (1.13), it can easily be seen that ¢(n) can be 
expressed as a generalized hypergeometric function (4.2) with 


¢¥(n) = Lin (1) 








(1.27) Dyo (p 

= 2 k=0 F7 

œo T(k+i)" 

= 2 k=0 FRF 

=) o kt L) 

= ¢(n, 1) 

= ¢(n) 
The value Ç” (0) =; Fo (1 |1) is singular and so must be calculated with the reflection equation 
(1.16) to get Lif (—1) =1 Fo( 1 |-1) = —4 = ¢(0) which agrees with the integral form of 
C(t)Vt #1 

oo sin(s arctan(s s?)72 . 

(1.28) Cli = (3 mt 2 fo es ds) ae Lio(—1) = -4 





2. NUMBER THEORY, CONTINUED FRACTIONS, AND FRACTAL STRINGS 


2.1. Fractal Strings and Dynamical Zeta Functions. A fractal string is defined as a nonempty 
open subset of the real line Q C R which can be expressed as a disjoint union of open intervals I; 
being the connected components of Q. [25, 3.1][30][23][22][11][24] 


(2.1) OE Sah 


The length of the j-th interval J; will be denoted @;. It will be assumed that Q is standard, meaning 
that its length is finite, and that £; is a nonnegative monotically decreasing sequence. 


lUa = VFL (E;)% < œd > 0 
Gh... >GG 0 





(2.2) 





where Jd > 0 means there is at least one value of d for which the statement is true. It can be the 
case that £; = 0 for some j in which case £; is a finite sequence. The sequence of lengths of the 
components of the fractal string is denoted by 

(2.3) £L = {4} 

The boundary of Q in R will be denoted by OQ = K C Q which will also denote the boundary of £. 
Any totally disconnected bounded perfect subset K C R, or generally, any compact subset K C R, 
can be represented as a string of finite length |Q|,. A subset K of a topological space Q is said to 
be perfect if it is closed and each of its points is a limit point. Since here Q is a metric space and 
K C Q is closed, the Cantor-Bendixon lemma states that there exists a perfect set P C K such 
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that K — P is a most countable. [35, 2.2 Ex17] As such, Q can be defined as the complenent of 
K in its closed convex hull, that is, Q = Q(K) is the smallest compact interval [a,b] containing K. 
The connected components of the bounded open set Q = (a, b)\K are the intervals I; of the fractal 
string £ associated with K. 





2.2. The Gauss Map h(x). Let Q, = (0,1)\OQn, where OQ, = {+00,0, 4 : n € Z} is the set of 
discontinous boundary points of the Gauss map h(a) € QaYx Z OQ, also known as the Gauss 
function or Gauss transformation, which maps unit intervals onto unit intervals and by iteration 
gives the continued fraction expansion of a real number 


Ma) =- [4] 
(2.4) = ne 


where |x] is the floor function, the greatest integer< x and {x} = x — |x] is the fractional part of 
x. [13, 2.1,3.9.1,9.1,9.3,9.7.1][42, A.1.7] Clearly h(x) is also defined outside of Q 


+ z>l 
(2.5) h(a)=4 h(x) -l<a<l 
+ +1 a4<-l1 
since 
1 0 xr>l1 
(2.6) B =4 || -l<e<l 
r -1 zr<-1 


As can be seen in Figure 1, h(x) is discontinuous at a countably infinite set of points of Lebesgue 
measure zero on its boundary 00», 





1 
(2.7) fu: lim h(x) # lim n(x) =O), = {400.0 =: n€ z} 
ry xr—>yt n 
The left and right limits of h(x) when x approaches an element on the boundary ôN, is given by 
limz-.aq- h(x) =0 
(2.8) lim, .aqt h(a) =1 


2.2.1. The Frobenius-Perron Transfer Operator. The Frobenius-Perron transfer operator [42, 2.4.4] 
[31, 2.3.3] [13, 1.3.1,8.2] [40, 1.8.2.4] of a unit interval mapping f(y) describes how a probablility 
density p(y) transforms under the action of the map. 


(2.9) [Usol(x) = f ôl- f(y))ely)dy 
where ô is the Dirac delta function and @ is the Heaviside step function. 
fo(a)dx = 0(x) 
(2.10) _ fO «<0 
A(z) ~~) 1 «20 


The function f(y) is the map being iterated and p(y) is some density on which the transfer operator 
U acts. Essentially, iteration of the map transforms points to points and iteration of the transfer 
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01 02 03 04 05 06 07 08 09 1.0 


— (hy K=O.) 


FIGURE 1. The Gauss Map 
operator maps point densities to point densities. The Gauss-Kuzmin-Wirsing(GKW) operator is 
obtained by applying the transfer operator to the Guass map. [44, 2] [50] [46] 


(2.11) De) =D, ae 


n=1 (n+x)? 





By changing the variables and order of integration in ne 37) an operator equation for ¢(s) is obtained. 





¢(s) awe [Una* "dx 
= s oTa- oie 
(2.12) =- sheli (x - (y u Ly* |) ye dydz 
= 2 1 GE) r dz 





An operator similiar to a is 


(2.13) (Sola) =D%2,0(2) - (ats) 
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The action of [U;,p](x) on the identity function x > x is given by 


ae 
[Cna\(z) = ra GaP 
z Da Gane 


WO (@+1) 
(2.14) = 


z 1 x£x+1 x+1 zxz+1 
re e+2 r74+2 2442 


(2+1)? 














where Y) (x) is the polygamma function (4.8). The area under the curve of [Unz] (£) is 





PO(r+1 1 
(2.15) Jo Wnal(a)dx =f. = Lape ; 
The identity action of [S;,p](a) is 
[Srz](z) = Xni E 
21 n=ln n+x 
(210) =7+U(¢+1) 


where y is Euler’s constant 


y = limno ea i — In(n) 
= lims ¢(s) — sI 
= lims41¢(s) + fE A(a)as~tda 
= lims» 4-8 J h(x)x® tdr + [P h(x)z®7tdz 
=1— J h(a)da 
~ 0.57721566490153286060651209 





(2.17) 


and the area under its curve is given by 
(2.18) SolSnz](a)de = f y+ (e+ de =1-7 
2.2.2. Piecewsise Integration of h(a). The Guass map h(x) € Qa is piecewise monotone [40, 2.1] 


between the points of OQ), and thus partitions the unit interval infinite covering set of decreasing 
open intervals seperated by Npa. [13, 5.7.1] Let In be an infinite set of open intervals 


(2.19) In =4 (sist) O<n<c 
(0,0) =0 n=co 
It is easy to see that 
(2.20) OQ, U Np = (0, 1] = Up=: Ty, 
[0,oo] = n=0 +n 





l-axn 1 1 
h (x) 2 T n+ cab si A 
(2.21) a 0 otherwise 


= (EP) (0 (SH) -9 (=) 
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where 6(t) is the Heaviside step function (2.10). We can reassemble all of the {h,,(a)}°2, to recover 
h(x) 











(2.22) U 
= De (At) (0 (258+) ae) 
where only one of the hp (x) is nonzero for each x. By setting n = |+| in (2.21) we get 
1-a| 4 | 1 1 
he) =4 A aS ST] 
(2.23) 0 otherwise 








= (S32) ARE) -°RP)) 
Define the partitioned integral operator |P f(a); a] (n) by 


(2.24) [Prea = PIE =P) = S7 fade = fy, fade 


n+l ty 


where by convention we have 


[Pf(z)](0) =lim,o+ fy, f(z)de = fy f(x)da 


oe (Pf(2)] (cc) =littns0 fi, fede =f? fade =0 
Thus 
(2.26) fo f@)dx =D PIE) = Cea Sr, f@)de 
i ieae SE PiOh@:.= >of t@ade 
Each interval J, has the length 
lIn = [P1] (n) 

(2.27) ao a 

~h p n+l 

= n(n+1) 


The elements n(n + 1) are known as the oblong numbers [15, A002378]. It is seen, together with 
(2.16), that 


In| = fg da 


a 1 
(2.28) = asi n(n+1) 


Í 
L 
à, 
= 


The piecewise integral operator [Pf(x); x] (n) can be used to calculate the area under the curve 
of h(x) which is also equal to the area under the curve of [S;,2](x). Let the length of the n-th 
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component hn(x) be denoted by 
hn = [Ph(x); 2] (n) 


(2.29) = f; ha(£)dz 


= fl hn(a)de 
— In(n+1)n4ln(n+1)—In(n)n—In(n)-1 
~~ n+l 








Regarding h(x) as a fractal string La = {h,(x)}2, its length |La] is given by 
Cn] = fo b(e)de 


7 So In(n+1)n+n(n+1)—In(n)n—In(n)— 1 





(2.30) nt] 

= fo (Snx](x)da 

= fy y+ U(a + 1)dx 
If n = 0 in (2.20) we get the interval Ip = (}, 4) = (1,00) and 
(2.31) tho= fp h(a)da = fi Eda = œ 


but if we choose a finite cutoff then 





(2.32) 1 = ialo) 
and 

œ ln(y) E 1 
(2.33) J y” dy E (n = 1)? 
thus 


” h(x)da oo oo fY tdg 
hax lore) Si h(a)d dy => fi zd dy 





n=2 1 y™ n=2 1 Le 
fore) co In 
E n=2 y Q) dy 
(2.34) see 
= ¢(2) 
re 
= 6 


2.2.3. The Mellin Transform. The Mellin transform [36, II.10.8][3, 3.6] is defined as 
(2.35) MEB f(x) = fe f(x)a®'dax 


where the usual definition of the Mellin transform is MLL) (x). Somewhat incredibly, by taking 


the Mellin transformation of h(x) over the unit interval, we get an analytic continuation of ¢(s) 
which is convergent when s is not equal to a negative integer, 0, or 1. When s is a negative integer 
or 0 the limit or analytic continuation must be taken since the series is formally divergent at these 
points, and of course the series s = 1 diverges. [45] [44] [43] 


My, h(t) = [Ph(x)a*~*; a] (n) 
(2.36) sje (2-|4]) a*-1d2 


n(n+1)78§+s(n+1)78§—n178 
s(s—1) 
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C(s) =2,-sM)h(2) 
=+- s fo h(x)as~*dax 
(2.37) =} -s fa (4-[4]) ote 





E co Ts 

~ ser" - Cp ee M;",5h(x) 

Sree oo n(n+1)78§+s(n+1)78—n178 
= 31 Sonan s(s—1) 














The term -#4 changes to -$4 = = — (—1) by subtracting the residue [47, 10.41][48, 6.1] of 








co 1s—1 
(2.38) ME h(x) = hle)z ae = f ay eee 
Io 1 


pepe ee 
s—1 





at the singular point s = 1, which happens to coincide with ee 





Res ( ite x de: 1) = Res (-4:1) 


x 


2.3. The Harmonic Sawtooth w(x). Define the harmonic sawtooth map w(x) € Qp,\OQp, which 
shares the same domain and boundary as the Gauss map h(x) to which it is similiar, and also 
has the property that its Mellin transform is the (appropriately scaled) zeta function. The n-th 
component w,,(z) is defined over the n-th interval In 


n+1 
0 otherwise 


=n(en+2—1) (o (2e) o (===) ) 
and by the substitution n > |4| we have 
w(x) = Dina Waha) ( ( ) ) 
(2.41) = 5 n(zn+ x -— 1) (6 (22+ o (==) 
=e beira) 


Unlike h(x) which is nonzero outside of |x| > 1, the (harmonic) sawtooth map has w(x) = 0VY|æ| > 1. 
The length of each component of w(x) is 


II 


= zi 1 
1) { n(an +a —1) Sman 


(2.40) 








(2.42) = fı, w(x)dz 





(2.43) = Leger fUn 
= Lin=1 atn 


NIK M 


TWO NEW ZETA CONSTANTS 


"| 


0.1 02 03 04 05 06 0.7 0.8 


—= mon 


FIGURE 2. The Harmonic Sawtooth 


The infinite set of Mellin transforms of w, (zx) 


Mie wl) = Mis wn (a) 
= [Pu(2)a*?; 2] (n) 
(2.44) = f™_ n(an +2 —1)2°~*dx 
. n+l 


0.9 


gzn—1 





= fy n(wn +x 1) (0 (=) -0 ( 
n(n+1)~8§+s(n+1)78—n1~ 
s(s—1) 





n 


1.0 


)) x®—lda 
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are summed to get an expression for {¢(s) : R(s) g No- } 


G(s) = s$ Jo LE] (e lz] +2- 1) 2 "de 

















=) n 18 sH Mir, w (x ) 
(2.45) =J s Se ce n(an + z — 1)as~1dax 
s+1 n(n+1)7$+s(n+1)7 -nts 
-Erns =) ) 


l-s =s 


ME n(n+1) $—n ~ $+sn 
= Xei s—l 





2.4. The Prime Numbers. Let P = {2,3,5,7,11, 13,17, 19, 23, 29,...} denote the set of prime 
numbers and N; = {1,2,...} and No = {0,1,2,...},N = {...,—2,—1,0,1,2,...} be the set of 
positive, non-negative, and signed integers. 

2.4.1. The Prime Counting Function: n(x). The prime counting function m(x) counts the number 
of primes less than a given number. It can written as 


(2.46) n(x) =>rPeP 4 


pcx 
which is essentially a step function which increases by 1 for each prime. [9, 15.11] 
2.4.2. von Mangoldt and Chebyshev’s Functions: A(x), 0(x), (x). Chebyshev’s function of the first 
kind @(x) is the sum of the logarithm of all primes < x 
(x) = ES m(pr) 
=In (sr) pr) 


where pẹ € P is the k-th prime. [7, 4.4] The generalization of m(x) is the Chebyshev function of 
the second kind 


(2.47) 


ve) = Eia h n(n) 
= llega x) | O(x E) 
(2.48) = ln(lcm(1, 2, sot: læ])) 
= Dn A(n) 
n(l—x~? =0 x? 
=g- ma-z) — ln(27) — ue $ ZVI) #0 
where the first sum ranges over the primes p € P and positive integers r and the sum over p is 


von Mangoldt’s formula where p ranges over the non-trivial roots of ¢(s) in increasing order. The 
function lem(....) is the least common multiple, and A(x) is the von Mangoldt function. 


= x 
Neo ae In(p) Lo pipe P,k € Ni} 
(2.49) 0 otherwise 


= Icm(1,2,...,7) 
=In (ES, ) 


A(s) is related to ¢(s) b 











d oo 
66 (5) A(n) 
2.50 — = VR(s) > 1 
(2.50) Stay eae VRE) 
Chebyshev proved that m(x), 6(a),andw(a) have the same scaled asymptotic limit. 
(2.51) i i ie A 
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[18, 1.3] [17, L4] [7, 4.3-4.4&3.1-3.2] [8] [9, 15.11]. Note that [9] incorrectly defines (x) as 
In(ged(...)). 


3. ANALYTIC CONTINUATION 


3.1. Continuation of ,,,F,, Near Unit Argument. The continuation formula for Gauss’s hy- 
pergeometric function 2Fı near unit argument is well known 





T'(a1)P' (a2) a, a2 _ œ (—1)”(1—z)” T(aitn)P(a2+n)0(s1—n) 
grt ( by z) = >D 





(3.1) T (bı) n=0 n! T'(a1+s1)P(a2+s1) 
s love) —1)"(1—z)” T(aitsi+n)l(ae+s1+n)C(—si-—n 
an ee aa ee, cD h Z a O ERA = 


where sı = bı — a, — a2 is the balance (4.3) of 2Fı which must not be equal to an integer, that 
is, 2F} cannot be s;-balanced. A function is said to be k-balanced only when k is an integer. When 
R(s1) > 0 the value at z = 1 is finite and given by the Gaussian summation formula 








ay ag 
2F1 b 
1 Ea T(bı—aı—a2) 
(3.2) T) Mayr ey 
SL 





T'(ai+s1)(a2+s1) 


3 ' It] = ¢*(1) = œ is 0-balanced and of course 
equal to the divergent harmonic series so the continuation formula does not apply. However, Biihring 
and Srivastava [6][5] generalized this relation to all ,41F, by expanding (3.2) as a series then 


interchanging the order of summations to derive a recurrence with respect to n 


It is obvious that lim;_,, Lif (t) = limy412Fi 








F Q1, .-.,ün+1 it _ T(bn)T(bn-1) 
3 3) n+l n by,..-, bn ~~ Pang)! (bn +bn-1-4n41) 
: oe (bn —An41)m(bn—1—Gn—1)m F Q1,--+,4n |t 
m=0 (bn +bn-1—-an+1)mm! n ay, bi, rao) bn—2; Dast + bn — An+1 +m 


which is valid V{R(a;) >0:1 <i < n+1}. The m-th term of the summand in (3.3) is contiguous 
(4.1.2) to the (m — 1)-th and (m + 1)-th terms and thus a linear relationship can always be found 
between neighboring terms. 





3.2. The Continuation of Li’ (t) and¢"(n) via Contiguous Functions. There are 4 functions 
contiguous (4.1.2) to Lif (t), only 3 of them are unique, and just 1 of them is interesting. The 
functions are obtained by shifting one of the numerator parameters a; + 1 or shifting one of the 
denominator parameters b; + 1. Shifting any of the a parameters or any of the b parameters will 
suffice since they are all equal and pF} is invariant with respect to the ordering of parameters. Let 
G@ and & denote the parameter vector €, where one element is shifted up or down by 1. 








GF =G-1,c+1l=6...,¢c+1 
3.4 Be ates oy 
ee) Ca = Enote n PS e Gem 7 


n—l1 
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For example, 43 = 4,4,5. Two of the four functions contiguous to Lif (t) are identical 
Hya Iny aE: 

(3.5) Badly 7 jt] t =n41 Fa a jt] t= Li; (t) 


Shifting any a; up is equivalent to shifting any b; down, both operations take Lif (t) to Lič; (t). 
Shifting any a; down results in the identity function since it puts a 0 in the numerator. 


(3.6) n+1Fn ( a t) t=t 
Thus, the only interesting function continguous to Lif (t) is obtained by shifting one of the denom- 
inator parameters up. Let this function be denoted by Li? F) 
1 = 
7 ly (2vt) — zh (2t) n=0 
(3.7) LIEH (t) =n41 Fa ( ntl t) =; $-—ł-1 n=1 


J+ 
an 





t 
CDe (1- BO DDH Lilt) n>2 


where I„(x) is a modified Bessel function of the first kind [34, 65] [10, 6.9.1] 


gh x 
(3.8) In(x) = Tia oft ( n+l =) 


Before applying (3.3), the notation will be simplified by extending (3.4) so that repeated shifts can 
be written more easily 


rea = a-1,C +j = Caer yQ CHj 
~x 


(3.9) e — 


5 : $ 
= Cn-1,6 7] 5 C,...,CGC— I 
SS 


n—-1 


where clearly Gt = @t! and & = @!. The goal is to extend Li¥ T! (t) to all Li” *™(t) by repeated 
application of Gree Applying (3.3) to (1.13) gives the continuation of Lif (t) > Lif, (t)Vn > 1 by 


setting @1,..n41 = In4qi and bi...n = > which results in 








(3 10) T tZ m=0 (m+1)(m+2) Vn >2 
& oF Jimy |t 
=t} m=0 (m-F1)(m-+2) Vn > 2 
since 





Plany) (bn +bn-1—any1) Ose = ty 2 
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and 


(bn — an4 1)m(bn 1— an 1m S (Da (Lm = 2 
(3.12) (bn + bn-1 — Gn41)m! = (2+2—1)mm! (m+1)(m+2) 








The denominator parameters 2,1} in (3.10) are simply 


DPN? sD 8 =2,...,2,3 
(3.13) n—1 25 +m ’ 94) +m 


n—2 


The numbers (m + 1)(m + 2) are known as the oblong numbers, [15, 4002378]. By simply setting 
t = 1 in (3.10) we get the continuation from ¢*(n) — F (n + 1)Vn > 1 





In 
F laii ee) ee Peas 
(3.14) ¢ (n) =n41 Fz ( Z, ) = ier (mt) (m3) Vn > 2 


The justification in saying that Lif (t) and ¢® (n) are continued to Lif ,(t) and ¢*(n+1) comes from 
the fact that the first term in the summand of the continuation (3.10) from Li; (t) + Lif (t) is 


; : : T, ' . : I 
contiguous to Lif; (t), that is, nFn—1 sei I) is contiguous to Lif_,(t) = nFn-1 ( = t); 
n—1 n—1 
The continuation formula (3.10) gives interesting answers for n = 0 and n = 1 which suggest an 
alternative to “the analytic continuation” of Ç (t) which is different from the usual —t— )>°_,(2n+ 


1—2-* 4n=0 
1)~*. We have 








oF, m+3 ) 
F = fore) 
oO 25s, \ a 
= ye — (Im 41(2)m+Im+41(2)—Im (2) T (m+3) 
mM (m+1)(m-+2) 
= Io(2) -1 
(3.15) x CE E 

o m+3 


CAL)! = Fog 86 | eee 





= x e(T(m+3)—-T(m+2,1)m—2r (m+2,1)) 
T m=0 (m+1)(m+2) 
= Fi(1) — 7 


= 1.3179021514544 


where Ei(x) is the exponential integral [10, 6.9.2] 


k 


Ei( ) E In(a~') ln(x) ae x 

Ke y 2 2 k=1 RT(kH1) 

(3.16) In(a—') In(x) 
=) 2 2 














1 1 
z2 Fo 2 2 |x 
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and T (a, z) is the incomplete Gamma function 


J a 
A ( a+1 |-2) 


a 





(3.17) T (a, z) =T (2) 


So we have the “hypergeometrically continued” values Ç (0) = Io(2) — 1 and ¢¥ (1) = Ei(1) — 7 
whereas the “real” values are ¢ (0) = — 4 and ¢ (1) = oo. In terms of reciprocal probability we have 


CF (0) = 78.15% 


(3.18) CE (1)! = 75.88% 


3.2.1. Lif (t) > Lif (t) and CF (t) > CF (t). The continuation ÇF (t) from n = 1 > 2 via (3.14) is 





straightforward 
1 11 
F 
E 3+) 
$ 1 1 
63 on 3+m 
a Amb (m+1)(m+2) 
(3.19) 


So, Nemt3 Met) 
Pas eae k=0 —T(m+k+3) 
> m=0 (m+1)(m+2) 








= fore) 1 (m+2) 
= 2 m=0 (m41) (m+2) (m+1) 
Fr 0 (m+1)? 


2 








6 


The continuation of Li? (t) to Lif (t) is a bit more complicated 











: 1 1 1 
Liz (t) = Fo ( 5 9 ) 
yid f 
(3.20) . ak 
= m=0 (m+1)(m+2) 
= at T2 (m, t) 
then rə(m, t) is given by 
pD (=1)P FIr (m+2)(Y(m-n+1)—4(m+2)) (1) e% (m+2);n 
rə(m, t) = neo nes nied) 
(3.21) (m+ 1)ePOrF2) tm+ 
(t-1)" + "77>" In(1—t) 





m+1 


so Lis (t) is equal to 
(3.22) 





ae m —1)"*1P(m+2)(U(m—n+1)—V(m+2 y( 1) Mehl) gh Aga, 
Lif (t)= y Nro = l - fp SS (t = Lyre 2 In(1 = t) 
2 or (m + 1)evm+2)em+1 m+1 
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where ¢)(m) = In(Icm(1, 2,3,...,m)) is Chebyshev’s function of the 2nd kind (2.48) and U(m) is 
the digamma function 
d del («) 


=— In(I'(x)) = The) 


(3.23) V(x) = 
3.2.2. ¢ (2) + ¢¥(3). The continuation from F (2) to Ç? (3) via (3.14) is carried out like so 
1 1 1 1 
F = 


(3.24) ee Ne 2 34m 
7 Emi (m+1)(m+2) 





T Den r3(m) 


Each term in the summand r3(m) has the form ao} +q3(m) where of course ¢(2 m and q3(m) 


) = 
is a rational function of m which follows a 3rd order linear recurrence equation[28, 8.2] given by 





qg(m) =q (m+ 1) (m? +8 m? + 21 m + 18) 
(3.25) q3 (m + 2) (—2 m? — 20 m? — 67 m — 75) 
34 


























q3 (m + 8) (m3 + 12 m? + 48 m + 64) 
—1 m=0 
(3.26) q3(m) = m=1 
08 m=2 
The solution to which is given by 
(42) — 
(3.27) gs(m) = SOO) 
so the summand r3(m) is 
(24m 
(3.28) r3(m) = nor + q3(m) = eee 
Thus (3.32) is also equal to 
co Ww (m+2 
(3.29) GE) = Ero EA 
Thus 
— WY (m+2) 
r3(m) — ati * 
— ¢(2,m+2) 
Zz m+1 
== T? m-1 1 
~~ 16 k=l kZ 
— Lk) Ga? 
g m+1 
(3.30) „[1 m+2 m+2 
ane m+3 m+3 
= (m+1)(m+2)? 
a 1 1 1 
e 2 3+m 





(m+1)(m+2) 
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The first 10 terms of {q3(m) :m=0...9} are 


5 49 205 5269 5369 266681 1077749 9778141 1968329 
8’ 108’ 576’ 18000’ 21600’ 1234800’ 5644800’ 57153600’ 12700800 


The denominator of (3.31) appears to be [15, A119936], the least common multiple of denominators 
of the rows of a certain triangle of rationals and the numerators are [15, A007406], the numerator 
of Xp ¢ from (4.9) which, according to a theorem Wolstenholme, p divides numer(q3(p — 1)) 
where p € P is prime. [12] [4] [1] 





(3.31) |—1, 


3.2.3. CE (3) + Ç (4). The continuation from F (3) to ÇF (4) via (3.14) is given by 


CE ooo a ) 


1111 
4F3 








(3.32) ae 2 2 34m 
= 2um=0 (m+1)(m+2) 
= J m=0 74(™) 
The summand r4(m) has the form 
(3.33) ra(m) =a(t, m) — b(t, m) — FeO te 4 CeO 


where a(t, m) is an (m + 1)-th degree polynomial and b(t, m) is a (m + 2)-th degree polynomial(the 
determination of which is left to an excercise for the reader or the topic of another article, but is 
readily obtained with the help of Maple[27]), and H (n) is the n-th Harmonic number 


H(n) = Xiz a 


=W(n+1)+y7 

(3.34) = J k1 eae 
2, 1 1 n+l 
ntg ? 2 n+2 


The polynomial b(t, m) vanishes when t = 1. An interesting set of formulas for ¢(4) is 


oS 2) (n4d)4 
ea) = O92, SRO) 


es (1) (m42) 
(3.35) yoo | YO (n+1) 42 Dio rt 
* 











n=1 2n(n+1) 


4. APPENDIX 


4.1. The Generalized Hypergeometric Function :, Fj. The Pochhammer symbol is defined according 
to 


T(n+k 
(4.1) (ne = TP 


The generalized hypergeometric function [39][48, 4.1] is defined as an infinite sum of quotients of 
finite products of Pochhammer symbols 


G1,...,4 o eT las) 
(4.2) oF, ( oe b ) = } k0 BT Gye 


-0q 
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The function ,Fj is said to be k-balanced [5] if the sum of the denominator parameters bı ...b, 
minus the sum of the numerator parameters a1 ...@)+1 is an integer. 


q Pp 
(4.3) k = bal(pFy) = X bn- X an 
n=1 n=1 


The value k is the characteristic exponent of the hypergeometric differential equation at unit argu- 
ment which is equal to the maximum root of the corresponding indicial equation and so determines 
the behaviour of the function near this point. A 1-balanced function is said to be Saalschiitzian. 
[41, 2.1.1] 


4.1.1. The Differential Equation and Convergence. The function ,F, converges when 


p<q V|t| 4 00 
(4.4) p=q+1 Vit} <1 

{p=q+1:ballpF,) > 1} vj=1 

p>qtl Vt =0 





where bal(,F,) = S<4_, bn — 5? _, an is the balance of the parameters (4.3). The differential 
přq n=1 


n=1°n 


equation solved by „F4 is of order max(p,q+1) 


q p 
(4.5) O: [[@+4 —1)-t][@ +a) f(t) =0 
j=l i=1 
where f(t) =p Fy ( n oF; es t) and @ = tå is the differential operator. When p = q + 1 
loses 
(4.5) has the form 


ae 


d 2 d 
(4.6) ao f(t) + af) P n - bn) F(t) = 0 


[48, 4.2][21, Ch3][34, 44-46] [41, 2.1.2] 


4.1.2. Contiguous Functions and Linear Relations. Any two hypergeometric functions 
pF4la...,b...;2) 
and „F4 (c...,d...;z) are said to be contiguous if all p + q pairs of parameters 


(a1,c1), BE. (ap, €p), (b1, d1), ee (bq, dq) 


are equal except for one pair which differs only by 1. There are 2p + q linearly independent 
relations between the (2p+ 2q) functions contiguous to ,F,(a...,6...; 2) where the relations are linear 
functions of z and polynomial functions of the parameters a., b... When any {a; = a,;:i 4 j} 
or {b; = bj :i Aj} in pF(a...,6..;z) there will fewer unique contiguous functions than if all the 
parameters were unique since the hypergeometric function is invariant with respect to the ordering 
of parameters. [41, 2.2.1] [34, 48] [39] [10, 4.3] [49] [33] 


4.2. Other Special Functions. 
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4.2.1. Polygamma Y™ (x). The polygamma function is the n-th derivative of the digamma (3.23) 
function 
(4.7) uw (r) = we) 

dx” 
and is defined as an infinite sum, a Hurwitz zeta function (1.4), or a hypergeometric function when 
x is positive integer 


co 1 1 
vz) = ( ae ae ma) A Same 
co n!(—1)” 
k=0 a nol 
1 12 
x—1 
n!( Ier: Tn++1 


=f e) n=0 


(-1)"*4*nld(2,n+1) n21 

or as a finite sum when zx is a positive integer and n = 1 [10, 1.16] 
(4.9) vO) =% -Pih 
4.3. Notation. 

Zo  {...,—2,—1,0,1,2,...} 

N o {1,2,3} 
(4.10) INS haan 

No  {0,1,2,3,....} 

No- {...,—3, —2, 1,0} 
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